The influence of polaron size on the conductivity of poly-DNA 
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The velocity of polaron migration in the long poly-DNA chain (~40 base pairs) in an applied 
electric field has been studied within a polaron model. We found that the polaron velocity strongly 
depends on the polaron size. A small polaron shows a slow propagation and strong tolerance to 
the electric field, while a large polaron is much faster and less stable with increasing electric field. 
Moreover, the conductance of the DNA molecule within the polaron model is found to be sensi- 
tive to structural disorders in the DNA geometry, but that dependence diminishes with increasing 
temperature. 



The application of the biological systems in develop- 
ing nanoelectronic devices has been recognized as one of 
the most intriguing and promising techniques in recent 
years. This is due to the molecular recognition and self- 
assembly properties that allows them to perform repara- 
tion of the damaged structures by the non-invasive tech- 
nology and self- incorporation of molecular blocks into the 
well-structured systems. The DNA molecule is one such 
biological system. The interest on DNA, while originally 
derived as it is the source of genetic information, has 
been significantly elevated due to the discovery of DNA 
conductance 0, 0, S II H 1 UMM- 

Understanding the mechanism of charge migration in 
a complex molecule such as DNA is a difficult problem. 
Just as in the field of condensed matter, in simple DNA 
sequences, if the donor and the acceptor are separated by 
a single potential barrier, the charge transfer has been ad- 
equately explained by the competition of quantum tun- 
neling and incoherent hopping 2|, Hi that has also been 
experimentally confirmed |lfj.lll|. However, the observed 
conductance of long DNA molecule with both poly- and 
mixed sequences, varies from insulator to a metallic be- 
havior 0,111, [EH- Further, the results for the decrease of 
conductance with decreasing temperature has been some- 
what contradictory 0, 0] ■ For the description of trans- 
port phenomena in polymer chains, a polaron model has 
been known to be very successful [12|, LUa] ■ For DNA, this 
model has shown recently to provide promising results 
as well, especially for explaining the temperature depen- 
dence of the DNA conductance 0, d, 11] • Moreover, addi- 
tional advantage of the polaron model, in comparison to 
the popular models such as the tight-binding approach 
or the system of kinetic equations, is that it includes the 
interaction of the migrating charge with the DNA lattice. 

In this paper, we have used the polaron model for sim- 
ulation of the dynamics of hole propagation in the dry 
poly-DNA molecules. We have found that the most prob- 
able reasons for variation of the DNA conductance in dif- 
ferent experiments 0, H, H, @] can be due to, (i) different 
external and internal conditions that can change forma- 
tions of large polarons to smaller ones, thereby decreas- 
ing the conductance. Therefore, the poly(dG)-poly(dC) 
and poly(dA)-poly(dT) chains due to different parame- 
ters for charge propagation should show different con- 
ductance behavior, as indicated in the experiment 0; 



(ii) Even a small disorder in the DNA structure can trap 
the polaron, again resulting in a decrease of the conduc- 
tance. However, a disorder-free DNA structure shows a 
constant velocity of polaron propagation in a constant 
electric field, which gradually increases with increasing 
electric field. 

We have simulated the polaron propagation through 
DNA using the Peyrard-Bishop-Holstein (PBH) model, 
which combines a quantum-mechanical treatment of 
charge motion with a classical treatment of the lattice 
distortion dynamics. Application of this model, which 
describes the charge tunneling from one DNA base site 
to another, is however limited to low temperatures. To 
avoid the temperature restriction we have investigated 
the poly-DNA molecule. In a poly-chain, each strand 
contains only one base type. Therefore, due to the large 
difference of ionization potential between the purine and 
pyrimidine bases 1.0 eV [14[) belonging to the oppo- 
site strands, the main mechanism of charge migration 
at low and at room temperature is longitudinal one- 
dimensional tunneling along a single strand containing 
purine bases. 

An electric field is applied to the system using the 
method of Ref. |13| . specifically, the charge transfer in- 
tegral V is multiplied by a complex exponential phase 
factor that includes the electromagnetic vector poten- 
tial A. For a uniform constant electric field E, we take 
A(t) = —cEt, thereby satisfying E = — dA/cdt. Here, 
c is the speed of light, and we have taken A=0 at time 
i=0. The uniform electric field is included via its uniform 
vector potential, rather than with its non-uniform scalar 
potential , in order that periodic boundary conditions 
can be utilized. 

With the modifications for the electric field, the cou- 
pled system of non-linear equations based on the PBH 
model is 
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FIG. 1: Polaron dynamics: the charge density \^n\ and the 
lattice displacement y n propagating through poly-DNA for 
£=0.016 mV/A(x=0.6 eV/A, V=0.1 eV, m=300 amu). 



where ^ n is the probability amplitude for the charge on 
the n-th base pair, V is the nearest-neighbor transfer in- 
tegral between the base pairs, x is the charge- vibrational 
coupling constant, e„ is the on-site energy for base site n, 
y n is the amount by which the n-th base site is displaced 
from its equilibrium position, to is the polaron mass on 
the single site, 7 is the friction parameter, Vm(2m) is the 
Morse potential, and W(y n , y n -i) is the interaction of 
neighboring stacked base-pairs, e is the electronic charge, 
and d is the interbase spacing (d=3.4A). The expressions 
and parameters for Vm(?m) and W(y n ,yi- n ) are taken 
from Ref. 0], and throughout this paper we take 7 = 1 
ps . In poly-DNA, x and V are constant with respect to 
the lattice site, and in the absence of disorder, we choose 
e n to be at the origin. 

In our simulations, a polaron is initially created by 
solving the system of nonlinear equations ([T]) and @ 
in a stationary situation (i.e., for d 2 y n /dt 2 = dy n /dt = 
and ihd^ n /dt = constant). Initial estimates for y n and 
\f„ are chosen to be nonzero only for ten consecutive base 
pairs. Due to the uniformity and stability of the system, 
the solution is largely insensitive to the initial estimates. 
The resulting initial solution can be seen in Fig. Q] for 
t = ps. The lattice displacement y n and the charge 
density \^ n \ 2 coincide and have the same overall shape. 
Together they comprise a polaron that, for the param- 
eters chosen here, is spread over approximately thirteen 
lattice sites but has a definite peak in the center. In 
the polaron's vicinity, the displacement y n , and there- 
fore the distortion energy XDn, are negative, creating a 
quantum well in the otherwise uniform energy profile of 
the poly-DNA. The wave function ^„ is consequently lo- 
calized within this well and the polaron is stable against 
small perturbations. Moreover, if no electric field is ap- 
plied (E=0), the polaron remains completely stationary 
because the initial solution corresponds to a stable equi- 
librium. 

In an electric field the positively charged polaron moves 



in the direction of the field (Fig. [T]) for a constant electric 
field of E=0.016 mV/A. Clearly, as the polaron migrates 
both the charge density and lattice displacement move in 
unison, i.e., the wave function and its formed quantum 
well travel together. Because of the stabilizing effect of 
the lattice distortion, the polaron retains the same overall 
shape as it migrates; however, our calculations indicate 
that in larger electric fields, the polaron becomes slightly 
more localized over time. 

In the following, we study polaron shape and prop- 
agation in poly-DNA under the influence of an electric 
field. We focus specifically on the effect of three physi- 
cal parameters, x, V, and to. The value of the charge- 
vibrational coupling constant x determines the decrease 
in on-site energy in the charged-state geometry, and has 
been estimated theoretically to be in the range of 0.3 - 
1.5 eV fl5| . Its value predominantly depends on the na- 
ture of the state geometry and its extension, which can 
be influenced by the structural parameters of DNA and 
the solvent environment as well. The charge in DNA can 
be spread in two direction: in parallel to the propaga- 
tion pathway - longitudinal direction and perpendicular 
- transverse direction. The spreading of the charge in the 
longitudinal direction significantly decreases x UM- For 
the transverse case, a charge can occupy a single purine 
base that suggests a larger value of x than does a charge 
that is partially delocalized over a base pair [l5| . The de- 
gree of charge delocalization also governs the value of the 
polaron mass m. If the charge is localized on a single base 
then to «150 amu, and this value increases with exten- 
sion of the geometry of the state in the transverse direc- 
tion. The value of the charge transfer integral V depends 
strongly on the configuration of the DNA geometry, and 
can be influenced by the solvent environment as well [l7| . 
Theoretical estimations place the value of V in the range 
of 0.05 - 0.3 eV MM, 

while experiments indicate the 
value of V to be <0.01 eV 0, 111]. Throughout this 
paper, we use the values x=0.6 eV/A, V—0.1 eV, and 
to=300 amu as reference points, and consider the effects 
of deviating from these values as we study polaron shape 
and polaron response to applied electric fields. 

The uniform system used in our simulations is ideal for 
studying the effect of model parameters on polaron shape 
since, as mentioned above, the polaron remains station- 
ary in the absence of an electric field. Figure [DJa) shows 
the effect of changing the coupling constant x under these 
circumstances. As is evident in the figure, increasing x 
results in greater polaron localization: the polaron oc- 
cupies fewer lattice sites, with a correspondingly larger 
charge density and lattice distortion at the center of the 
polaron. Conversely, decreasing x has the opposite ef- 
fect. Figure [2^b) shows the shape of the polaron at vary- 
ing values of the charge transfer integral V. The effect is 
opposite to that observed when adjusting x : increasing 
(decreasing) V causes a decrease (increase) in polaron lo- 
calization. This is expected since the larger the charge 
transfer integral, the more the polaron will spread out to 
neighboring sites. In accordance with Eq. @, the value 
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FIG. 2: Polaron shape in poly-DNA for varying values of (a) 
X and (b) V. In (a), V=0.1 eV, m=300 amu, and x is given 
in units of eV/A. In (b), x = 0-6 eV/A, m=300 amu, and V is 
given in units of eV. 



of m has no effect on the polaron's shape: in the station- 
ary state the time derivative my dy n /dt in Eq. ([2]) equals 
zero. 

After the polaron is initially created, a constant uni- 
form electric held E is applied and we study the result- 
ing polaron motion in time. Simulations show that there 
exists a maximum electric held E max that the polaron 
can tolerate; for E > E max , the lattice displacement 
and charge density become unsynchronized and irregular. 
The value of E max depends on the system's parameters 
and in all cases studied, the more delocalized the po- 
laron the smaller the value of E max . As such, increasing 
X or decreasing V leads to a polaron capable of tolerating 
larger helds, whereas E max is mostly independent of m. 
A highly localized polaron causes a large localized lat- 
tice distortion, which effectively acts as a quantum well. 
By decreasing the energy in its vicinity, a more local- 
ized polaron is therefore more stable, and consequently, 
it is reasonable that the polaron would be able to remain 
well-formed in the presence of stronger electric fields. 

The application of an electric field effects polaron mi- 
gration (Fig. [T]). Moreover, provided that E < E max , 
the polaron remains well-formed and the motion is con- 
tinuous. Further, in the presence of a field the polaron 
moves at a constant velocity. This result can be explained 
as follows: (i) our system represents poly-DNA under pe- 
riodic boundary conditions, so the system parameters do 
not vary with location; (ii) the applied electric field is uni- 
form and constant; (iii) the polaron maintains a constant 
shape as it travels; and (iv) the friction term rwy dy n /dt 
prevents the charge from continually accelerating. These 
factors ensure that after a brief period of acceleration 
from the initial stationary state, the polaron's velocity 
does not vary. 

The velocity does depend, however, on the system pa- 
rameters and on the applied electric field strength. Fig- 
ures GHa), G3h), and [3]Jc) show how the velocity of the 
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FIG. 3: Propagation velocity of a polaron through poly-DNA 
as a function of applied electric field for varying values of (a) 
X, (b) V, and (c) m. Points corresponding to E > E max 
are absent as they do not represent accurate values. Unless 
specified otherwise in the figure, X = 0-6 eV/A, V=0.1 eV, and 
m=300 amu. 



polaron varies with E for different values of \, V, an d 
m respectively. These graphs all display a similar behav- 
ior, but the magnitudes and slopes in the plots are very 
sensitive to the values of Xj V, an d m. In particular, an 
increase in the velocity — and therefore conductance — 
can be achieved by decreasing the coupling constant x> 
by increasing the charge transfer integral V , or by de- 
creasing the polaron mass m. Analysis of both \ and V 
suggest that an increase in polaron localization produces 
a decrease in velocity. It should be noted, however, that 
this correlation need not be a general mathematical rule; 
for example, by changing m it is possible to change the 
polaron's velocity without affecting its localization at all. 
Other model parameters, such as the friction constant 7, 
also affect the velocity, but do not change the qualitative 
behavior observed in Fig. [3] 

Structural disorder in the poly-DNA can drastically 
alter the polaron propagation. An irregularity at a base 
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FIG. 4: Effect of an irregularity: the charge density \^ n \ 2 and 
the site energy xVn + £n propagating through poly-DNA with 
a disorder-induced well. Here £20=0.1 eV and e n =0 for n^20. 
At t=0 ps, X2/20=0 eV; at t=100 ps, xy20=0.029 eV. 

site k can produce a quantum barrier or well, which we 
represent with a non-zero value of efc. If efc > (a bar- 
rier), the polaron is usually unable to traverse site k for 
any E < E max , instead coming to a stop just before the 
irregularity. Only for very low barriers (efc «0.01 eV) in 
conjunction with the polaron size and a large electric field 
(E raO.l mV/A) can the large polaron migrate through 
the barrier. The situation for quantum wells (efc <0) is 
very different. For very shallow wells (e^ >-0.05 eV), the 
polaron completely enters the well for any E, although a 
large E can dislodge the polaron if the well is extremely 
shallow. For deeper wells, the polaron sometimes splits 
into two, with one polaron halting prior to the well and 
second entering the well (Fig. [3}. In these cases, the ini- 



tial polaron stops before the irregularity, but then part of 
it tunnels from the polaron- induced well (xHn) into the 
disorder- induced well (efc). It should be emphasized that 
the above analysis applies only to poly-DNA. In less uni- 
form situations, the system is generally unstable, and the 
polaron often tunnels through barriers to reach distant 
wells 0. 

In conclusion, the velocity of polaron propagation and 
hence the DNA conductance is mostly determined by the 
polaron size in transverse and longitudinal directions. A 
decrease of x induced by partial derealization of the 
charge from the purine base to the whole base pair [3| 
provides the extension of polaron size in the longitudi- 
nal direction and significantly increases the DNA con- 
ductance. For example, a decrease of \ by 0.2 eV can 
increase the conductance by ~ 6 - 10 times depending on 
the electric field value. The polaron can be destroyed at 
high electric fields, which then changes the mechanism of 
charge transfer in the DNA molecule and therefore, will 
cause a discontinuity at the conductance characteristics. 
A large polaron has higher velocity but is less tolerant to 
the electric field. The structural disorders, which form a 
barrier < 0.05eV (well) on the polaron pathway, depend- 
ing on energetic conditions and electric field magnitude, 
can stop (trap) the polaron and cause a fast decrease 
of DNA conductance resulting in an insulating behav- 
ior. The influence of disorder on the DNA conductance 
decreases with increasing temperature. 
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